Atom-laser coherence and its control via feedback 
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We present a quantum-mechanical treatment of the coherence properties of a single-mode atom 
laser. Specifically, we focus on the quantum phase noise of the atomic field as expressed by the 
first-order coherence function, for which we derive analytical expressions in various regimes. The 
decay of this function is characterized by the coherence time, or its reciprocal, the linewidth. A 
crucial contributor to the linewidth is the collisional interaction of the atoms. We find four distinct 
regimes for the linewidth with increasing interaction strength. These range from the standard 
laser linewidth, through quadratic and linear regimes, to another constant regime due to quantum 
revivals of the coherence function. The laser output is only coherent (Bose degenerate) up to the 
linear regime. However, we show that application of a quantum nondemolition measurement and 
feedback scheme will increase, by many orders of magnitude, the range of interaction strengths for 
which it remains coherent. 

PACS numbers: 03.75.Fi, 42.50.Lc, 03.75.Be 
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I. INTRODUCTION 

The invention of the laser in the late 1950s [[j] created 
the field of quantum optics and continues to lead to an 
enormous range of scientific and technological applica- 
tions. It is expected that the realization of "atom lasers" 
will similarly revolutionize the field of atom optics [Q. 
Atom optics is the study of atoms where their wavelike 
nature becomes important, suggesting an analogy with 
photons ||. An atom laser is therefore defined as a de- 
vice that produces a continuous beam of intense, highly 
directional, and coherent matter waves 0, in analogy 
with the light produced by an optical laser . The ideal 
atom laser beam is a single frequency (i.e., monochro- 
matic) de Broglie wave with well-defined intensity and 
phase. 

The first experimental achievements of the Bose- 
Einstein condensation (BEC) of gaseous atoms || was 
followed immediately by several independent ideas for 
creating an atom laser |Q. Since then there have been 
experimental advances in the coherent release of pulses 
H and quasicontinuous beams || of atoms from BECs, 
as well as further theoretical proposals pQ| . In the 
experimental configurations to date, the laser mode is 
the ground state of a trapped BEC, which is pumped 
by evaporative cooling of uncondensed atoms, and the 
out-coupling (separated in time from the pumping) is 
achieved by either Raman or radio-frequency (rf) transi- 
tions to an untrapped state. Although these experimen- 
tal accomplishments do not include simultaneous pump- 
ing and output coupling, they do represent major steps 
towards achieving an operating atom laser. 

Recent experimental (ll], |l2|, [l3| and theoretical || [l4|, 
|l5| , [l6[ [It], Q [l9| studies have focused on the funda- 
mental coherence properties of BECs and atom lasers. 
Atoms (unlike photons) interact with each other, pro- 
ducing strong nonlinearities that affect the coherence of 
the trapped condensate and thus also the out-coupled 



laser field. For a single-mode condensate the dominant 
effect of atomic collisions is to turn number fluctuations 
into fluctuations in the energy and hence fluctuations in 
the frequency, thus causing increased phase uncertainty. 
Collisional interactions therefore lead to a significant de- 
crease in the atom-laser coherence time, and a corre- 
sponding increase in the linewidth, especially in the case 
of BECs formed by evaporative cooling. However, we 
have previously shown that a continuous, quantum non- 
demolition (QND) feedback scheme can effectively cancel 
the linewidth broadening due to such collisions [^0| . 

To study the coherence properties of an atom laser 
one can either focus on classical or quantum noise in the 
atomic field. Sources of classical noise may be technical, 
such as fluctuations in the trapping potential, finite tem- 
peratures, and specific trap geometrie s |12|| , or dynami- 
cal, such as three-body recombination ]19{. The study of 
these effects is usually based on mean-field laser models 
described by Gross-Pitaevskii (GP)-type equations 
Quantum noise is an intrinsic part of the atomic system 
as a consequence of the uncertainty relations and is the 
limiting contribution to coherence. The study of this 
requires a fully quantum-mechanical approach, of which 
the most common is based on the quantum optical mas- 
ter equation [ p2| . The complexities of either approach, 
which individually require approximations to facilitate 
theoretical analysis, indicate that simultaneous analysis 
would not be easy p3| . 

In this paper we present a fully quantum mechani- 
cal treatment of the coherence of a single-mode atom 
laser and its control via the QND feedback scheme pro- 
posed in Ref. (^0|. Specifically, we study the properties 
of the first-order coherence function (i.e. phase coher- 
ence) , which allows us to derive the coherence time (and 
hence laser linewidth) as well as the power spectrum of 
the laser output. Section II summarizes a set of require- 
ments for the coherence of an atom laser first detailed in 
Ref. B . Section III presents our mathematical model for 
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the atom laser and shows the resulting linewidth for in- 
creasing atomic interaction strength. Numerical methods 
are presented in Sec. Ill B, while the analytical results 
are discussed in Sec. Ill C. When the collisional nonlin- 
earity is very strong, the coherence function undergoes 
a quantum collapse and revival sequence (detailed in III 
D). This leads to an interesting regime in the power spec- 
trum (detailed in III E), which has not been considered 
before. Section IV details the effects of feedback based on 
a physically reasonable QND measurement of the conden- 
sate number, to include all regimes of the nonlinearity. 
Section V concludes. 



II. REQUIREMENTS FOR COHERENCE 

The coherence of an atom-laser beam can be defined 
analogously to that of an optical laser beam Q. The 
fundamental assumption is that the laser output is well 
approximated by a highly directional classical wave of 
fixed intensity and phase, which is also ideally restricted 
to a single transverse mode. The output should also be a 
stationary process, i.e., its statistics should be indepen- 
dent of time. To be coherent, the laser output should 
then additionally have (1) a relatively small spread of 
longitudinal spatial frequencies (i.e. be monochromatic); 
(2) a relatively stable intensity (i.e. be approximately 
second-order coherent); and (3) a relatively stable phase 
(i.e. have a relatively slow decay of first-order coherence) . 

The first condition, monochromaticity, follows from 
the requirement that the laser output approximates a 
classical wave. It can, therefore, be expressed in terms 
of the characteristic coherence length of the wave, Z co h = 
(<5fc) _1 , i.e. the reciprocal of the spatial frequency spread. 
Condition (I) becomes 

?coh > A = 27r/fc, (2.1) 

i.e., that the coherence length be much greater than the 
mean atomic wavelength [H. In terms of the spectral 
width or linewidth, £ = 8lo, the monochromaticity re- 
quirement simply becomes £ <C where the mean fre- 
quency Q is defined by the kinetic energy of the atoms, 
i.e., Q = fik 2 /2m. Condition (I) also guarantees that the 
dispersion of an atomic beam will be negligible over the 
coherence length [|J . 

To explain the second and third conditions we require 
a many-body description of the output beam; see Ref. Q 
for an in-depth discussion. Basically, the output field of 
a laser can be represented by the localized field anni- 
hilation operator b{t), which approximately satisfies the 
(5-function commutation relation 

[b(t),tf(t')] = 6(t-t'), (2.2) 

at a given point in the output. Thus, I{t) = (t)b{t) 
can be interpreted as the approximate atom-flux opera- 
tor. The fundamental assumption of a laser is that the 
output should be well approximated by a classical wave 



of fixed intensity and phase. This is therefore represented 
mathematically by 

b{t) » P{t) = 0e- iQ \ (2.3) 

where (3 is a complex number and the trivial time depen- 
dence emphasizes that the laser output should be sta- 
tionary. The atomic field is not exactly a classical wave 
because there are fluctuations in the field amplitude due 
to classical and quantum sources of noise. These will need 
to be small or somehow cancelled to maximize coherence. 

It can be argued that there are no mean fields in both 
atom and quantum optics p4| , and as such the descrip- 
tion of a laser as a coherent state is a "convenient fiction" . 
Despite this, mean- field theories (e.g., based on GP equa- 
tions) are very successful in describing the properties of 
lasers, as is the use of initial coherent states for solving 
and visualizing master equations (as done in this paper). 
A vanishing mean field means that 

(&(*)) = 0. (2.4) 

In the case of optical lasers the output might, in principle, 
be in a state with well-defined amplitude and phase. But 
since we do not know the absolute phase, an average over 
all possible phases gives Eq. ( |2.4| ). For atom lasers, only 
bilinear combinations of the atom field are observable 
and similarly no Hamiltonian is linear in the atom field 
(atoms cannot be created out of nothing). Thus a mean 
field amplitude is physically impossible and the absolute 
phase is unobservable, again giving Eq. ( |2.4| ). 

Since the mean field of a laser is zero, we cannot require 
{b{t)) = (3. However, for approximating &(£) by Pit), we 
can require that the mean intensity be given by 

(J(i)) = <&t(t)6(t)) = |/3| 2 , (2.5) 

and also that the fluctuations in intensity should be 
small in some sense. This requirement is quantified us- 
ing Glauber's normalized second-order coherence func- 
tion (for a stationary system) p5| : 

g^{r) = {:I{t + r)Iit):)/(Iit)) 2 , (2.6) 

where the (: :) denotes normal ordering. 

For a field that is second-order coherent, i.e. g^(r) = 
1, there is no correlation between the arrival times of 
bosons at a detector and their distribution is Poissonian. 
Specifically the probability for detecting a boson in the 
interval it + t, t + r + dt) given one detected at time t 
is g^ 2 '(r) (lit)) dt. For the intensity fluctuations to be 
small we therefore require that jj] 

k/ 2 >(r)-l|«l, (2.7) 

i.e. the laser output should be approximately second- 
order coherent: condition (II). 

Assuming that condition (II) is met, the intensity of 
the laser beam will be relatively stable and the only sig- 
nificant variation in the output field will be due to phase 
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fluctuations. A useful measure of the phase fluctuations 
is the stationary first-order coherence function (25): 

GW(r) = (6t(t + T )6(t)> ) (2.8) 

or its normalized form: 

g^{T)=G^{T)/(b\t)b{t)). (2.9) 

Unlike the field b(t) itself, the bilinear combinations 
above are measurable even for an atom field. G^(t) 



is simply the mean intensity (2.5) when r — and as r 
increases, it decreases towards zero as the phase becomes 
decorrelated from its initial value at t. 

So, the phase of the field might be undefined because 
it varies in time, i.e., the first-order coherence decays. 
Although we cannot expect the laser to be approximately 
first-order coherent for all time [i.e., <7^(t) ~ 1], we can 
require the decay of (r) to be slow in some sense. The 
characteristic time for this decay is simply the coherence 
time, which can be defined as 
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9 {1) (r)\dT. 



This is the time over which the phase of the field is rela- 
tively constant. 

But even if the phase is constant, it might also be 
undefined due to a large intrinsic quantum uncertainty 
given by AfiAn > 1/2 |2{|. Since typically An < n, 
the quantum phase uncertainty will be large if the mean 
number n is small. For the phase to be well defined we 
therefore need the field to have a large intensity, i.e., n ^> 
1, over the time that the phase of the field is constant, 
i.e., for times T <C T co h- The number of bosons in the 
output field for a given duration T is h — (I{t)) T. Thus, 
for a well-defined phase we require (I(t)) r co h 3> 1- In 
terms of the first-order coherence function this translates 
to | 

(/(*)) / \g<»{T)\dT = [ \G^(r)\dr » 1, (2.11) 



which quantifies the requirement that the decay of first- 
order coherence be relatively slow: condition (III). 

This condition for coherence is equivalent to the 
requirement that the output field be highly Bose- 
degenerate Q, which is rarely considered for optical 
lasers because it is so easily satisfied. It requires that 
the output atom flux, i.e., (I(t)), be much larger than 
the linewidth 5u>. Since the atom-laser linewidth is the 
reciprocal of the coherence time, i.e., 



1= Suj = l/r< 



Coh; 



(2.12) 



we find that condition (I) requires £ <C u> and condition 
(III) requires 



£ « </(*)> 



(2.13) 



For single-mode optical lasers typically (I(t)) > u> so con- 
dition (III) is always satisfied. On the other hand, the 



collisional interactions in atom lasers causes significant 
linewidth broadening and so satisfying (III) is not guar- 
anteed. However, in Sec. IV we show that this broaden- 
ing can be effectively cancelled by a QND measurement 
and feedback scheme. 

The linewidth, or spectral width, of a field is usu- 
ally defined as the full width at half-maximum height 
(FWHM) of the output power spectrum. In general, the 
power spectrum is given by the Fourier transform of the 
first-order coherence function |E2|: 



P(w) 



G (1) (r)e" 



r dr. 



(2.14) 



This is defined so that J™ P{uj)duj = (I(t)) and there- 
fore can be interpreted as the steady-state flux per unit 
frequency p7f . 

If the first-order coherence function has the form 
G^'(t) oc exp(— 7t) then the laser output will have a 
Lorentzian power spectrum. In this case the FWMH is 



exactly equal to the linewidth as defined in Eq. ( 2.12 ), i.e. 
the reciprocal of the coherence time. On the other hand, 
(/.±UJ j£ G^(t) oc exp(— 7 2 t 2 ) then the laser has a Gaussian 
power spectrum with a FWHM that is now only approx- 
imately equal to the linewidth. 

In any case, if the first-order coherence function has the 



form GW(t) = |G (1) (t)I exp(iwr), then condition (III) 
for coherence can be restated in ter ms o f the maxi mum 
spectral intensity P(u>). From Eqs. (2.14) and (2.10) and 



the above assumption (Appendix A shows that this is a 
good approximation for the atom laser), we find 



P(co) 



\G^(r)\dT = 4(I(t))T< 



coh ■ 



(2.15) 



Thus, regardless of the resultant shape of the output 
power spectrum, condition (III) becomes 



P(Q) > 1. 



(2.16) 



Note that the central frequency ui will be shifted by any 
laser dynamics which cause a rotation of the mean phase 
of the laser field. 

The remaining sections of this paper present a study 
of the quantum phase dynamics of the atom laser and 
thus the coherence properties of its output. Specifically 
we study the first-order coherence function as a measure 
of phase fluctuations, which (in the absence of intensity 
fluctuations) will be the limiting factor to the coherence 
time of an atom laser. This function, as indicated in this 
section, is also intrinsically related to the laser output 
characteristics of linewidth and power spectrum. 



III. ATOM-LASER LINEWIDTH 
A. Atom-laser dynamics 

The atom-laser model consists of a source of atoms ir- 
reversibly coupled to a laser mode, which is supported in 
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a trap that allows an output beam to form. A broadband 
reservoir acting both as a pump and a sink is also coupled 
to the source modes. The laser mode and source can be 
modelled by the ground and excited states of a trapped 
boson field, representing the condensed and uncondensed 
atoms, respectively 0J. Gain can be achieved, for exam- 
ple, by evaporative cooling of the uncondensed atoms. 
Out-coupling from the laser mode can be accomplished, 
for example, by coherently driving condensed atoms into 
an untrapped electronic state. This model can be simply 
described by a quantum optical master equation for the 
laser mode alone ^, ^| , obtained by adiabatically elimi- 
nating the source modes and tracing over the continuum 
of output modes. The system is thus characterized by 
a wavefunction <&(]") and annihilation operator a for the 
condensate mode. 

Far above threshold, the laser mode has Poissonian 
number statistics ]3(^, |3l|. In the absence of thermal 
or other excess noise, its dynamics are modelled by the 
completely positive master equation 0, B9J 



p = KfxT>[a/]A[(v] 1 p + kV[o\p = Cop, 



(3.1) 



where the superoperators T> and A are defined as usual 
for an arbitrary operator r: 

V[r]p = rpr 1 - A[r]p, A[r]p = |{r f r, p}. (3.2) 

That the master equation is of the Lindblad form follows 
from the identity 



V[a^}A[a^} 



where C is known as the Liouvillian for the total system 
evolution. 

As master equations, Eqs. ( |3.l| ) and (3.5) are de- 
rived using the Born-Markov approximation. More com- 
plicated mathematical models for the atom laser dy- 
namics may include a full multimode description with 
non-Markovian pumping and/or damping (see, e.g., 
Refs. 0, 18 1). However, only the single-mode master 
equation description employed in this paper allows a rel- 
atively straightfor ward analysis. Although a single-mode 
scheme, e.g., Eq. (3.5), ignores the many source modes 
and the continuum of output modes, it is the simplest 
physically reasonable model for an atom laser, in that it 
includes the essential mechanisms of gain, loss and self- 
interactions. 

There are also physical justifications for using Marko- 
vian theory. Reference JlTj states that the Markov ap- 
proximation is valid for output coupling rates satisfying 
T m , where T m is the output memory time, which 
typically ranges from 10~ 2 to 1 ms. This condition is 
thus easily satisfied for typical values of which range 



from 10 -2 to 10" 



Furthermore, as discussed in 



dqV[a} e- qaa ^ 2 ]. (3.3) 



Ref. |Ti|, the Born-Markov approximation is only valid 
for either weak output coupling or large atomic densi- 
ties. BECs formed by evaporative cooling have strong 
atom-atom interactions that correspond to large atomic 
densities. We are therefore just ified in making the Born- 
Markov approximation in Eq. (3.5) [but not necessarily 
in Eq. (3.1)]. In other words, we expect that strong non- 
linear interactions, rather than any non-Markovian dy- 
namics, will dominate the linewidth. 



The first term of Eq. (3.1) represents linear output 
coupling at rate k and the second term represents non- 
linear (saturated) pumping far above threshold, where 
p S> 1 is the stationary mean boson number. It is the 
decreasing difference between the gain and loss, as the 
laser is pumped above threshold, that gives rise to the 
gain- narrowed laser linewidth j p2| . The localized output 
mode operator b of the preceding section is then related 
to the laser mode via b = v + -y//ca, where v represents 
vacuum fluctuations fl32|] . Note also that we have cho- 
sen a reference potential energy for the system such that 
there is no Hamiltonian oc a^a in the master equation. 

To include the effects of atom-atom interactions in the 
laser mode we consider a simple s-wave scattering model 
for two-body collisions, which is valid for low tempera- 
tures and densities [[33]]. This is described by the Hamil- 
tonian: 



Hroii = hCataJaa, C 



|$(r)| 4 ^r, (3.4) 



where $(r) is the condensate wavefunction and a s is the 
s-wave scattering length. The total master equation for 
the laser mode including atomic interactions is then 

p = KpV[a i ]A[a ! ]~ 1 p + nV[a]p — iC[cr a* aa, p] = Cp, 

(3.5) 



B. Numerical calculation of linewidth 

As stated in the Sec. II, the coherence time of a laser, 
r co h, is roughly the time for the phase of the field to 
become uncorrelated with its initial value. As shown by 
Eq. ( 2.10D , it is dete rmi ned by the stationary first-order 
coherence function (2.9), in which the output operators 
b can be replaced by the laser mode operators a, since 
b = v + */Ra. For the evolution described by the master 
equation (|3.5|), the coherence function becomes 



9 {1 Ht) 



Tr[a f e £t ap s 



]/Tr[A 



.at a] 



(3.6) 



where p ss is the stationary solution to Eq. (3.5) given by 

mi): 



\re i0 )(re i6 \, (3.7) 



where r = y/Ji. Thus, the state of the laser (3/7) can 
be thought of either as a mixture of number states or 
equivalently a mixture of coherent states (although see 
p5| for a discussion of this). 

It is a very good approximation (see Appendix A) to 
assume that if g^\t) is not real then its complex nature 
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is simply of the form g^\t) = \g^' (t)\e iut where Q is 
the central frequency of the laser output. That is, we 
assume g^(t) is complex due to an effective detuning 
C a = — iQ\a*a, p] in the evolution. This type of evolution 
causes a rotation of the mean phase proportional to Co, 
whereas th e de cay of |<?^(t)| indicates phase diffusi on. 
Using Eq. ( |3.6| ), this allows the integral in Eq. (2.10) to 
be evaluated to give 



7"coh 



-Tr[a t (£-iw) 1 ap ss ]/2Tr[p s 



,a'a| 



(3.8) 



Equation (3 



can be evaluated numerically, for ex- 
ample using the Matlab quantum optics toolbox [ p6| . 
The first guess for Co is found from the approximation 



Im{Tr[a^£ap ss ]} 



wo, 



(3.9) 



which is exact for C = C a . Subsequent corrections 
are foun d by an iterative procedure. Substituting tok in 
Eq. ( |3.8| ) gives tu, which is then used to update our guess 
for Co via the expression u>k+i = tOk — Im(l/2rfc). Basi- 
cally, this scheme ensures that the calculated r co h has a 
vanishing imaginary component. This is ju stifie d in Ap- 
pendix B, where we also show that, if Eq. (3.8) is valid, 



then only one correction to coq is needed for an accurate 
determination of Co. 



C. Analytical calculation of linewidth 



Analytically, it is easier to use the fact that Eq. (3.6) is 
unchanged if p ss is replaced by the initial coherent state 
\re ){re \ for arbitrary 9 (say 8 = 0). We therefore have 

gM(t) = TV[otp(t)]/r, p(t) = e ct \r)(r\. (3.10) 



Using any suitable phase-space (a, a* 
this expression is then equivalent to 



gV>(t) = (a*(t))/{a*(p)), 



representation, 



(3.11) 



where (a*(0)) = (a(0)) = r and \a) is a coherent eigen- 
state of the laser field, i.e., a\a) = a\a). The state of the 
field at any time can thus be described by the probability 



distribution for a, or equivalently (since a 



') the 



intensity, n = \a\ , and phase, tp, distributions. 

The fluctuations in intensity are relatively small for a 
laser with /i > 1, i.e., Sn(t) ~ 0. Then, also assuming 
the number statistics are unchanged by the evolution, we 
have n(t) w n(t) = n(0), which gives 



9 {1) (t) 



,-itp(t)\ 



,-itp(t)\ 



(v^(0)e-^( )) (e-'V(°))' 



(3-12) 



Now the phase distribution at time t due to the laser 
evolution is given by ip(t) = ip(Q) + (p{t), i.e., the phase 
distribution of the initial coherent state plus the rela- 
tive phase difference <p{t) = arg[a(i)/a(0)]. Assuming 



this phase evolution is independent of the initial phase 
uncertainty, we have 



(3.13) 



where the second approximation assu mes Gaussian 
statistics. The coherence time, Eq. ( 2.10D , is thus found 
by evaluating the integral: 
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(3.14) 



In the Q- function representation [|22[ , a density opera- 
tor p has a corresponding Q function defined by 



Q(a,a*) = (a\p\a)/ir, 



(3.15) 



normalized such that J Q(a,a*)d 2 a = 1. The action 
of an operator on p thus has the corresponding mirror 
action of a differential operator acting on Q(a, a*). This 
is the most convenient representation for our laser model 
because of the identity (see Appendix C) 



E 

fc=i 



d_ 

dn 



Q(n,p), 



(3.16) 



which, since the higher order derivatives are negligible, 
can be truncated at k — 2. This representation also al- 
lows us to visu alize the dynamicsproduced by the master 
equation (3Ji) as shown i n F ig, ffl. 

The master equation ( |3.5[ ) thus turns into a Fokker- 
Planck equation (FPE) for Q(n 7 ip): 



St 



Q(z,t) 



d 2 



3,k 



^— [Aj(z)Q(z,t)], (3.17) 



where z = (n, <p) and the drift vector A and diffusion 
matrix B are given by 



A = 



K(p 
(3- 



-1 -n) 
2n)C 



B = 



2/t(/i + n) 
2nC 



2nC 
k/2u 



(3-18) 

To find equations of motion for the moments (zj) and 
(zjZk) the FPE needs to be converted to an Ornstein- 
Uhlcnbeck (OU) equation. For an OU process, the drift 
vector is linear in the variables (n, ip) and the diffusion 
matrix constant. Our drift vector A is already linear, 
but our B matrix is not constant. The simplest option 
is thus to replace all the amplitudes in B with their (Q- 
function) mean value, i.e., n — * p + 1. The equations of 
motion for the moments are then d(zj) jdt = (Aj) and 
d (zjZk) /dt= (zjA k ) + (z k Aj) + (Bjk + B kj )/2, where A 
and B are now OU parameters. 

We find that the number statistics are unchanged from 
that of the initial coherent state, which for the Q-function 
representation are fi = p, + 1 and V n = 2/i + 1. However, 
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FIG. 1: Typical contour plots of the Q function for fi — 15 
and C(= «x/4/u) = k,(2ty u)^ 1 ^ 2 . Solid ring: stationary laser 
state p ss , given by Eq. (3.7). Solid circle: initial coherent 
state \y/ fi + +T|; dashed contour: phase diffusion due 

to laser gain and loss at t ~ 0.8/k; dotted conto ur: state due 
to total evolution including collisions, i.e., Eq. (3.5), also at 
t ~ 0.8/k. 



the phase-related moments are altered, 
they are (for the relative phase <j>) 



4>{t) 



-2/j.Ct, 
8/iC 2 



e" Kt + Kt - 1) 



K 
K 



Using /x 3> 1, 



(3.19) 
+ ±t, (3.20) 

(3.21) 



As expected, there is a mean phase shift ( [3.19] ) due to the 
collisions, while the nonzero covariance (C n ^ = {n<j)) — 
(4>)) given by Eq. ( |3.21 ) explicitly shows the number- 



phase correlation produced by the collisions. The phase 



variance (3.20) contains two terms, where the second cor- 



responds to standard laser phase diffusion. The first term 
thus indicates the increase in phase fluctuations due to 
collisions. 

The effect of collisions on the phase fluctuations can 
be clearly seen in Fig. [l]. This figure shows single con- 
tour plots of the Q-function for the hypothetical coherent 
state | •»//!} and snapshots at a later time due to the evo- 
lution of the master equation ([T^). If we ignore collisions 
the effect of the laser evolution is simply phase diffusion. 
By including collisions we see two effects. Firs t, there is 
a rotation of the mean phase, due to Eq. ( 3.19| ), and sec- 
ond there is phase shearing. T his is due to the nonzero 
number-phase correlation, Eq. (3.21), indicating that if 
the inherent number fluctuations produce n > n the cor- 
responding phase 4> will be less than <fi, and vice versa. 
The initial co here nt state will approach the actual laser 
state p ss , Eq. (3.7), as t — > oo. 

Substituting Eq. (3.20) into the express ion fo r magni- 
tude of the first order coherence function (3.13) gives 



\9 {1) (t)\ 



l)/4/*—«i/4ji 
6 i 



(3.22) 



where we have introduced x = ^hC/k as a dimension- 
less parameter for the atomic interaction strength. This 
expression does not have a simple analytical solution. 



However, by inspection, there are two limits that can 
be solved analytically. If \ <C JJi we obtain 



2r, 



coh 



Mx 2 +iW^ dt = 4^/«(l + x 2 ). (3.23) 



For x yfu , on the other hand, the first exponential 



in Eq. (3.22) is dominant and then expanding e 



to 



second order we obtain 

/>oo 

2r coh a / e-* 2 * 2 * 2 / 8 ^ = ^/2^/k X - (3.24) 
Jo 

The resultant expression for the atom laser linewidth 
due to collisions is thus 



l= { ^(l + X 2 )/2/i for X < 
2kx/ v / 2tt^ for X > 



(3.25) 



Clearly, for y <C 1 , we obtain the standard laser linewidth 
&o = k/2^ |29[ These two expressions agree at 

X — y/8[i/ir, and they are an excellent fit to the nu- 
merical calculations of Eq. ( |3.8| ) , except at the boundary 
between the regimes. This is illustrated by the figure in 
our previous paper p0| , and also the extended version in 
this paper, F ig. p | appearing in Sec. Ill D. 

Equation ( |3.25[ ) represents the same physical dynam- 
ics as found in similar studies b y t he authors of Refs. 
and [^|. Zobay and Meystre jlS] present a three- mode 
atom-laser model, with the output mode adiabatically 
eliminated. Ignoring collisions between pump and laser 
modes, they obtain phase variances [Eqs. (21) and (22) 
of Ref. [15|1, which are similar and identical to the ex- 



ponents of Eqs. ( 3.23| ) and (3.24) respectively. See also 
the linewidth plotted in Fig. 3 of Ref. p5| . Gardiner 
and Zoller p^J studied a Bose-Einstein condensate in dy- 
namical equilibrium with thermal atoms. Our first-order 



coherence function, Eq. (3.22), has the same structure as 
the analogous expression Eq (184), derived in Ref. [§7|. 
The two regimes of Eq. ( 3.25 ) correspond to the charac- 
teristic time constants of Eqs. (187) and (186) of Ref. |37j] 
respectively. The second of these expressions, where the 
nonlinearity is dominant, is familiar as the inverse col- 
lapse time of an initial coherent state in the absence of 
pumping or damping |38, p9[. 



Since the output power spectrum is the Fourier trans- 
form of G^(t) [see Eq. (2.14)], the shape of the spec- 
trum is also determined by the form of V$ (t) . For the two 
regimes of Eq. (3.25), the laser output has Lorentzian and 
Gaussian power spectra respectively, as was also found 
in Ref. ||l5| . These spectra are illustrated by Fig. || in 
Sec. IV B. See Sec. Ill C for a more in-depth discussion 
of the atom laser power spectrum. 

The standard laser linewidth (in the absence of col- 
lisions) is simply given by re/2/x. For the preliminary 
atom laser experiments of Refs. J8|, |9|, the interaction 
strength C is always found to satisfy C»k//i and hence 
X ^> 1 [Q. Atom lasers, therefore, have a linewidth far 
above the standard limit. Furthermore, if x ~ M 3 ^ 2 the 
linewidth will be larger than the mean output flux Kfi. 
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In other words, the atomic collision strength does not 
have to become very large before the laser outp ut do cs 
not satisfy condition (III) for coherence [i.e., Eq. (2.13)]. 
It is thus of great interest to find methods for reducing 
the linewidth due to atomic interactions. One method is 
continuous QND measurement and feedback as shown in 
Sec. IV. 



D. Revivals of the coherence function 

In the preceding section, the atom-laser linewidth was 
calculated for atom-atom interactions ranging from weak 
(x 1) to strong (x ~> \fP)- However, exact numerical 
calculations based on Eq. ( |3.8| ) indicated that there is an 
upper bound to the linewidth (occurring for \ > p 2 ) that 
was not included in the previous analysis. It turns out 
that, in addition to linewidth broadening, the collisional 
interactions also lead to quantum revivals |38| of the first 
order coherence function. Although, note that in this 
very strong collisional regime, the output atomic beam 
cannot be considered a laser according to the definitions 
in Sec. II (since £ > Kfi for x M 3 ^ 2 )- 

To study the regime of revivals it is helpful to start 
by ignoring all other laser dynamics apart from the col- 
lisions. In this case, p = — iC[a^a^aa]p = CqP, and we 
can analytically solve for the periodic structure of g^\t). 
For an initial coherent state p(0) — \a)(a\, 



Using the number state representation, i.e. pit) = 
J2Pmn(t)\n)(m\ and \a) = exp(-|a| 2 /2) J2 a n \n}/s/rd, 
we find for the first-order coherence function 

9 «(t)=exp[-/i(l-e 2iC % (3.27) 

since \a\ 2 = p for the laser, and its magnitude is 

| 5 (1) (i)| = exp[-^(l - cos 2(7*)]. (3.28) 



The coherence function clearly has periodic revivals when 
t = rrm/C, m is an integer. 

Including the other laser dynamics, i.e. gain and loss, is 
not so straightforward. Since these terms, unlike those in 
Cc, are not functions of the number operator we cannot 
easily utilize the number state representation. The main 
effect, however, is si mply a decaying envelope applied to 
the revivals of Eq. (3.28), such that the strength of C 
compared to np will determine the number of significant 
revivals in the coherence function. This will be shown 
below. The revivals of the coherence function become 
significant when x ~ 47r/i 2 or C k kit p. This regime 
was determined by calculating t he e xact linewidth based 
on numerical solutions of Eq. (3.8) and corresponds to 
the interaction strength where the linewidth begins to 
approach a maximum. 

To determine the value of this maximum linewidth, we 
extend the work of Milburn and Holmes (40| , who model 
an anharmonic oscillator coupled to a zero-temperature 
heat bath, via two basic assumptions for including sat- 
urated gain. The master equation modelled by Milburn 
and Holmes is (using our notation) 

Pmh = -iC[(o t o) 2 ,/9] + KV[a]p, (3.29) 

which gives a first-order coherence function of the form 



gW(t)=Tr[<Jp(t)]/a; p(t) = e c ^\a)(a\. (3.26) g^ n (t) 



e lCt e- 



-Kt/2 



exp 



p(l - in/2C) 



(1 



e 2lCt e 



1 + k 2 /4C 2 

(3.30) 

To add saturated gain to this model we first assume 
that, far above threshold, the contribution to phase diffu- 
sion is equal for both gain and loss and so we replace 
k with 2k. Second, including gain will almost cancel the 
overall exponential decay e~ Kt l 2 of the coherence due to 
loss, resulting in the smaller term e ~ Kt / 4 i- 1 (since this will 
give the standard laser linewidth n/2p). These assump- 
tions give the following results for g^(t) and ^^(i)!: 



9 (1) (t) 
l.9 (1) WI 



e~ Kt/i ^ exp 



p(l 



1 + K 2 /C 2 



1 + K 2 /C 2 



We only expect this expression to be valid in the very 
strong interaction regime (C 3> up). Here revivals of 
|5^(t)| are signifi cant and k/C<1 (since p ^> 1). Also, 
as shown by Eq. ( 3.2S ), revivals occur at mt r = mir/C, 
so in the strong regime the envelope of the coherence 



JliCt „—1kV\ 



: cos2Ct - — e 

c 



i2Ct. 



(3.31) 
(3.32) 



I 

function is given by (for finite m 
|<? (1) (*)U = |<7 (1) (mt r )|~exp 

~ e- 2t/t «, t Q = 1/ up. 

Here the time tQ can be interpreted as the quantum dis- 
sipation time, i.e. how long the state would last if it was 



Ap 



+ 2np I mt r 



(3.33) 
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in a superposition of coherent states. Nonunitary effects, 
such as damping, cause a decay of the quantum coher- 
ence of these states at a rate cx Kfi [[l2| . This is relevant 
because the state produced halfway between revivals by 
nonlinear interactions such as Cq is in fact a superposi- 
tion of coherent states f43|| . 

The relationship between the quantum dissipation 
time and the revival time can be used to give an indi- 
cation of the number of significant revivals in the coher- 
ence function for a given interaction strength. That is, if 
t r ^> tQ then no revivals will be seen, but if t r <C tQ as 
for the above equation, the number of significant revivals 
is of order tQ/t r . Revivals begin to appear at t r ~ £q, 
which is at \ — 47r/i 2 or C ~ k-k\x as stated earlier. 

We are now in a position to determine coherenc e tim e 
and linewidth in the regime of revivals. From Eq. (2.10), 
the coherence time is simply half the area under the 
function \g (1) {t)\. In the very strong interaction regime 
(C 3> K/i), this area will be made up of many individual 
peaks w hich decrease in height due to the envelope given 
by Eq. ( 3.33|) . The first of these peaks (which actually 
starts at |g^(0)| = 1) will be the same as the coherence 
function for no revivals, and its area will be ^/2tFJ1/kx 
as given by Eq. ( p. 24 ). The subsequent peaks will have 
areas twice this area multiplied by the height of the en- 
velope at that time. Thus we have for the total area 




-2t/t Q 



^2 S (t ~ mt r )dt - - > 

m=0 J 



m— 



e -2mt T /t Q _ ±_ 
2 



(3.34) 



This expression can be evaluated by using t r <C tQ and 
by noting that for a geometric series, J2m rm = V(l — r )- 
The analytical expression for the linewidth in the regime 
of revivals is then 



4kv27t^ 



,3/2 



(3.35) 



This is in exact a gree ment with the numerical results 
obtained from Eq. ([D]), which are plotted Fig. |2|. 

In this figure we have plotted the approximate an- 
alytical expressions for the linewidth in the absence 



[Eqs. fl3.25| ) and ( p.35[ )] and presence [Eq. (|yj)] of feed- 
back (see Sec. IV B for details) as a function of the non- 
linearity x f° r M — 60. We have also included numeri- 
cal results (see Sec. Ill B for details) as a comparative 
test for the analytical work, which is valid for /x ^> 1. 
As can be seen, the agreement is very good even for an 
occupation number of only 60 thus confirming the 
accuracy of our analytical expressions for the linewidth. 
Without feedback we see four distinct regimes. There is 
the standard laser linewidth for x « 1, a quadratic de- 
pendence on x f° r 1 X ^ v^*' an< ^ a l mear regime for 
s/Ji <C % <C /i 2 . The latter two correspond to the regimes 



of Eq 



.25!) . Finally there is a constant regime given by 



Eq. (3.35|) for X ^ A* i which is due to the collapses and 




X = 4uC/k 10 

FIG. 2: Atom laser linewidth in units of k for rj = 1 and 
fi — 60, plotted with and without feedback using both ana- 
lytical (lines) and numerical (points) methods. The dotted 
line corresponds to £ = k/j,, i.e. interaction strengths with 
corresponding linewidths below this line satisfy the coherence 
condition of Bose degeneracy [condition (III)]. For the feed- 
back results, see Section IV B. 



Note that the ap pro ximation used in the numerical 



calculation of Eq. (^8|), i.e., g^(t) = \g (1) (t)\e""-, is 
no longer necessarily valid in the regime of revivals, 
as indicated by the multi-complex-exponential nature of 
Eq. flQgj ). Nevertheless, our numerical results are still 
correct because we have taken the mean atom number fi 
to be an integer. At revivals the approximation to (t) 
becomes 



\g^\mt r )\e lUomU = \g^ \e i2m7T ^ , 



(3.36) 



revivals of the coherence function. 



where we have only used the first guess for u>, since the 
iterative procedure will be inaccurate in this regime (see 
Appendix B). This expression clearly equals \g^\ for in- 
teger [i. Thus, the same numerical simulation can be 
used for all values of C as long as fj, is an integer 



E. Power spectrum 

As stated in Sec. II, condition (III) for coherence re- 
quires that the i ntegral of | G^ 1 ^ (i) | be much greater than 
unity [Eq. (|2.1l[)]. This was reinterpreted as requiring the 
linewidth to be much less than the output flux (Bose de- 
generacy), or equivalently requiring the maximum spec- 
tral intensity to be much greater than unity P{ui) 3> 1. 
Now the linewidth is only the FWHM of the power spec - 
trum if it is Lorentzian. As discussed after Eq. ( 3.25 ), 
this will only be the case in the weak-interaction regime 
X *C y 7 ^. As x ( or C) is increased the power spectrum 
becomes Gaussian, and as x enters the very strong in- 
teraction regime it will no longer have a simple structure 
at all. At these strong values of the nonlinearity we see 
quantum revivals of g^\t). 

In terms of normalized first-order coherence function 
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FIG. 3: Output power spectrum P(lo) where uj has units of 
k, plotted for /i = 15 and various values of C in the strong 
interaction regime. The dotted line corresponds to the value 
4K/^/^ max , i.e., the maximum height of the power spectrum in 
the regime of revivals. Note the log scale for uj. 



the power spectrum becomes 

r- OO 

P(uj) = K fi g^\t)e~ lu)t dt, (3.37) 



where we have recognized that (/) = (i^fr) = k (a^a) = 
K/i. From this equation it is clear that as long as g^(t) 
has a simple structure, i.e., no revivals, then the spectrum 
will have a simple (Lorenztian or Gaussian) lineshape for 
a given interaction strength C, with the intensity and 
width determined by how fast g^(t) decays. 

The maximum spectral intensity was defined in Sec. II 
by Eq. ( p 15h , which was based on the assumption that 



Since we are in a reference frame 



with zero mean frequency before including collisions, uj is 
the detuning frequency due to collisions which causes the 
mean phase shift of an initial coherent state (see Fig. [I]). 
At this frequency, the output power spectrum will, there- 
fore, have a maximum value given by 



P{uS) = 4K[lT coh . 



(3.38) 



As stated above, in the regime of revivals this approxi- 
mation will only be accurate for the first guess for uj, i.e., 
luq = 2\iC. 

The linewidth in the strong atomic interaction regime 
(C > K7r/i) was calculated in the preceding section to 
be AkV2tt^ 3/2 , i.e. Eq. ( |3.35| ). Thus, from the above 
equation, we expect the maximum spectral intensity to 
reach 



P(uj ) = 4Kfi/t max = lh 



(3.39) 



as the atomic interaction strength C is increased far 
above unfi. 

Figure [§] illustrates both the simple lineshape of the 
power spectrum when there are no revivals and the com- 
plication of the spectrum as C is increased. Since all 
the plotted spectra have P(u~i) < 1, the atom-laser out- 
put c learl y does not satisfy condition (III) for coherence, 
Eq. (2.16), in the strong interaction regime. The first 



spectrum in Fig. || is in the range ^/Ji < C < nirfj,, and 
thus although revivals are not seen, the output is still 
above the cutoff for Bose degeneracy. 

The remaining spectra illustrate the increasing effect of 
quantum revivals due to increasing interaction strength. 
The central peak of the power spectrum, which c an b e 
defined regardless of revivals as shown in Eq. (2.15), 
clearly approaches the predicted maximum of 1/ y/2nfi, 
i.e., Eq. (3.39), for C S> nnfi. For examples of the spec- 
trum in the weak interaction regimes of Eq. (3.25), see 
Sec. IV B. 



IV. REDUCING THE LINEWIDTH VIA 
FEEDBACK 



Section III A showed that the atomic interactions do 
not directly cause phase diffusion. Rather, they cause 
a shearing of the field in phase space, with higher am- 
plitude fields having higher energy and hence rotating 
faster. The resultant linewidth broadening is a known 
effect for optical lasers with a Kerr (x' 3 ') medium [ E5[ . 
The shearing of the field is manifest in th e fini te value ac- 
quired by the covariance Cn^it) in Eq. ( |3.2l| ). It means 
that information about the condensate number is also 
information about the condensate phase. Hence, we can 
expect that feedback based on atom number measure- 
ments could enable the phase dynamics to be controlled, 
and the linewidth reduced. 



A. QND feedback scheme 

QND atom number measurements can be performed 
on the condensate in situ via the homodyne detection of 
a far-detuned probe field Q ^ |8[ fl|] . This dispersive 
interaction causes a phase shift of the probe proportional 
to the number of atoms in the condensate. We consider a 
far-detuned probe laser beam of frequency uj p and cross- 
sectional area A that passes through the condensate |30| . 
Figure |] shows an experimental schematic for our QND 
measurement and feedback scheme. 

For a single atom the interaction with the optical probe 
field can be approximated by the Hamiltonian 



V 



4A \8AAI sa J PP 1 Pl 



(4.1) 



where £l p , A, 7 have their usual meaning and / sa t = 
2nhc^/\ 3 |53j|. Here p is the annihilation operator for 
the probe beam, normalized so that huj p p'p is the beam 
power. 

The effective interaction Hamiltonian for the whole 
condensate can thus be taken to be 



-Hint = h6(a*a - n)p^p, 



(4.2) 



where 6, defined in Eq. (4.1), is the phase shift of the 



probe field due to a single atom. Here we have also sub- 
tracted the mean phase shift such that the probe optical 
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condensate "anti-mean" lens 



probe 



B-field 



l * ; 8<|>ac8« 




FIG. 4: Experimental schematic. A far-detuned probe laser 
of amplitude g interacts with the condensate (i.e., the laser 
cavity) causing a phase shift in the probe proportional to 
the condensate number. The "anti-mean" lens subtracts the 
mean phase shift leaving the probe with a phase shift pro- 
portional to the number fluctuations Sn. The photocurrent 
Ihom(t) from the homodyne Y detection of this field is then 
used to modulate an external _B-field uniformly applied to the 
condensate. 



laser is a measure of the atom number fluctuations only 
(see Fig. |). 

The back action on the condensate due to this interac- 
tion can be evaluated using the techniques of Sec. Ill B of 
Ref. 0| . Assuming the input probe field is in a coherent 
state of amplitude g and mean power P, the evolution of 
the atomic system due the measurement is 



P = Q 2 T>[e 



-i${ar a— fj.)] 



MV[a f a]p, 



where the the measurement strength is given by 



M = q z 9 z = P6 z /hu) p . 



(4.3) 



(4.4) 



The approximation in Eq. ( |4.3| ) requires 9{a* a — fx) <C 
1, which for Poissonian number fluctuations is simply 

V/i0<i. 

The above result represents decoherence of the atom 
laser due to photon number fluctuations in the probe 
field, resulting in increased phase noise. In a recent theo- 
retical study by Dalvit and co-workers , it was shown 
that dispersive meas ure ments of BECs cause both phase 
diffusion [as in Eq. (O)] and atom losses. Nevertheless, 
they also show that phase diffusion dominates the de- 
coherence rate for large atom numbers, i.e., for /i » 1, 
and so the depletion co ntri bution can generally be ne- 
glected [[33|. Equation (<h3) is equivalent to the corre- 
sponding phase diffusion term in their work, since it can 
be shown that both M and the phase diffusion rate given 
by Eq. (14) of g§ reduce to ~ y*\ 5 I/hcAA 2 . 



The effect of the interaction (4.2) on the output probe 
field is to cause a phase shift proportional to the number 
fluctuations, a^a — p. The output field operator is given 
by H 



Poat = e v 



'pn 



' Pn 



iQ9{a)a-p), (4.5) 



where again the approximation requires ^fpZB <C 1. Ho- 
modyne detection of the Y quadrature of the output 
probe field will thus be a measure of the condensate num- 
ber fluctuations. The homodyne photocurrent operator 
is given by p2[ 



2\/M(a t a-^). (4.6) 




FIG. 5: Contour plots of the Q function for p, = 10 and 
C = k(2typ)~ 1 ^ 2 (ignoring any mean phase shifts). Black 
circle: initial coherent state | ^/Jt) {^/p\ ; black ring: stationary 
laser state p ss , Eq. (B.7). The other contours correspond to 
the evolution at time t ~ 1.5/ k due to successive terms in 
the master equation. Dashed contour: phase diffusion due 
to £0; dotted contour: including atomic collisions, i.e., £0 
and C; dot-dash contour: including QND measurement back 
action, i.e., £0 and C and M, and finally the solid contour 
corresponds to Eq. ([f ,9|). All contours are for the optimal 
feedback regime F — ^fr\M = C. 



In order to control the phase dynamics of the conden- 
sate, we wish to use this homodyne current to modulate 
its energy. This can be done, for example, by applying 
a uniform magnetic field or far-detuned light field across 
the whole condensate. In the ideal limit of instantaneous 
feedback, we model this by the Hamiltonian 



H & (t) = -ha f aFI hom (t)/VM, 



(4.7) 



where F is the feedback strength and Ihom(i) is the clas- 
sical photocurrent corresponding to the operator lj ut . 

The total evolution of the system including feedback is 
obtained by applying the Markovian theory of Ref. |32| . 
The master equation becomes 



p = CqP — iC[a^ a) aa, p] + MT>[a^ a]p 
+ iF\a} a) aa, p] H T>\a)a]p, 



(4.8) 



where we have allowed for a detection efficiency 77 |52| ] 
and dropped terms corresponding to a frequency shift. 

The terms in Eq. (4.8) describe, respectively, the stan- 
dard laser gain and loss (Co), the collisional interactions 
(C), the measurement back action (M), the feedback 
phase alteration (F), and the noise introduced by the 
feedback. As before we can visualize the effect of the 
measurement and feedback terms on the evolution of an 
arbitrary coherent state. This is illustrated by the Q- 
function contours in Fig. [5]. In this figure we have ig- 
nored the mean phase shift due to collisions to make the 
comparison clearer. 

To completely remove that unwanted nonlinearity, the 
obvious choice for the feedback strength is F = C . We 
also want to minimize the phase diffusion introduced by 
both the measurement and feedback. A weak measure- 
ment will give poor information about the atom number, 
with a high noise-to-signal ratio, which will increase the 
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noise due to feedback. On the other hand, if the measure- 
ment is too strong the measurement back action itself will 
dominate. This leads us to guess the optimal regime for 
both measurement and feedback to be F = ^JrjM = C, 
which simply leaves 

2C 

p = KuV\a^]A[a)}- 1 p + KVlalp + — V[a'a\p. (4.9) 

Vv 



B. Linewidth results 

Proceeding as before, we can find the exact effect of 
the general feedback scheme on the atom-laser linewidth. 
Neither the measurement nor the feedback affect the 
atom number statistics. The change in the phase statis- 
tics are reflected by the new Fokker-Planck equation for 
Q{n,<j)). The altered terms in the drift vector and dif- 
fusion matrix are A 2 = (3 — 2n)(C — F), B\ 2 — B 2 i = 
2n(C - F), and B 22 = n/2n + M + F 2 /rjAl. After lin- 
earizing, we again find the phase-related moments: 



C n< p(t) 



-2p{C - F)t, 
8/i(C-F) 2 , 



Kt- 1 



2p 

MC~F) 



F 2 \ 



1) 



(4.10) 

(4.11) 
(4.12) 



where again the approximations have used p 1. 

These equations clearly show that all the unwanted 
phase statistics are cancelled by choosing a feedback 
regime with F = C and furthermore the minimum phase 
variance is when M = C/^/rj. Specifically, both the 
mean phase shift and the correlation between number 
and phase fluctuations are removed, and the phase vari- 
ance is simply given by 



f k_ 2C 
V2M + 



t. 



(4.13) 



This is ex actly the phase variance from the master equa- 
tion (U). 



In this case, Eq. (3.14) has a simple analytical solution: 



'coh 



2// 



s/v) 



(4.14) 



where we have again used the dimensionless atomic inter- 
action strength \ — ^pC/n. Note that unlike Eq. (3.25), 
this linewidth is valid for all \- The derivation of 



Eq. (4.14) above is based on preselecting the feedback and 
measurement parameters (i.e., F — ^/rjM = C). On the 
other hand, Ref. |20) presented the analytical solution for 
the linewidth independent of the choice of these parame- 
ters and proceeded to find the minimum with respect to 
the feedback strength. The result [Eq. (24) of Ref. |o|] 
only differs from Eq. (4.14) by the term k(— l/Arj)/2p. 
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FIG. 6: Output power spectrum P{ui) for p — 15 vs u> in 
units of k, plotted with and without feedback for interaction 



strengths of \ 
■ A 4 



3 ' 2 and 



X 



p . The dotted spectrum 



(x = p 3 ^ 2 with feedback) has a maximum intensity of ~ 31. 



Fr om F ig, g (in Sec. Ill C) and the above result, 
Eq. (4.14), it is evident that our QND feedback scheme 
offers a linewidth much smaller than that without feed- 
back for most values of x- ln fact f° r \fji- ^5 X ^ 47r^ 2 
the reduction in linewidth due to feedback is a factor 
of \J%pj-K. Most importantly, the laser output is Bose 
degenerate [satisfies condition (III) for coherence], up to 
X ~ p 2 with feedback, as opposed to x ~ A* 3 ^ 2 m the 
absence of feedback. Thus, the atom laser with feed- 
back remains coherent for much stronger atomic non- 
linearities than without feedback. It is interesting that 

this corresponds to the "conditionally coherent" regime 
^3/2 < ^ < ^ 2 ag di scussec i i n R e f Q _ 

Since the nonlinearity C is effectively cancelled b y thi s 
feedback scheme, the output power spectrum, Eq. ( [3.37 ), 
will never have the compl icate d structure as shown in 
Sec. Ill C. Also, since Eq. (4.13) has a linear dependence 
on time (rather than higher powers) the first-order co- 
herence function decays exponentially and as such will 
produce a Lorentzian output power spectrum [see dis- 
cussion after Eq. (2.14)]. 

This lineshape is illustrated in Fig. ||, which also plots 
the corresponding spectra that would be produced by the 
laser (for the same values of x) without feedback. These 
latter sp ectra have a Gaussian lineshape a s disc ussed af- 
ter Eq. (|3.25| ). Also, as indicated by Eq. (|4.10|) , the ro- 
tation of the mean phase due to collision is cancelled 
by our feedback scheme, i.e., g^'(t) is no longer com- 
plex and hence Co — 0. The power spectrum for the atom 
laser including feedback will thus be centered around zero 
frequency regardless of the atomic interaction strength, 
which is confirmed in the figure. 

The spectra in Fig. [fj] are plotted at the two cutoff 
values for Bose degeneracy, which are x ~ A* 3 ^ 2 for no 
feedback and x ~ A* 2 when feedback is included. If the 
output is coherent, i.e., satisfies Bose degeneracy, then 
it clearly has a much narrower and thus more intense 
spectrum than if the output is not coherent. This figure 
thus illustrates the link between the three definitions of 



linewidth [discussed after Eq. ( 2.15 )] and its relation to 
coherence. That is, a narrow intense output spectral line 
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corresponds to a long phase coherence time, and more 
specifically, if the FWHM linewidth is much less than 
the value of the output flux, then the laser output will 
satisfy the conditions of coherence. 



C. Experimental realizability 

We will now briefly examine some issues of experimen- 
tal realizability. A question of interest is, how easy is 
it to obtain a QND measurement of sufficient strength 
to optimize the feedback? We have shown in Sec. IV B 
that the optimum feedback scheme requires N = C'/^/rj. 
From Ref. p4| ] we know that most current experiments 
work in the regime where the Thomas-Fermi approxima- 
tion can be made, allowing an analytical expression to be 
found for \ [Eq. (6.25) in Q] and C. Typical values are 
X ~ 10 3 , n ~ 10 6 and C ~ lCTV 1 . To determine 6 and 
hence M we use typical 87 Rb imaging parameters |f|, fl9| , 
which include A = 780nm, A ~ lCr n m 2 , 7 — 5Mhz, 
A ~ 2GHz, and J sat ~ 10W/m 2 . For these values, 
9 ~ 3.3 x 1CT 6 , and thus M ~ 4.2 x 1CT 4 x I. To obtain 
a measurement strength of the order of C, we therefore 
require a probe laser intensity of only ~ 30W/m 2 , which 
is quite reasonable. 

A related question is, how much of a problem is atom 
loss due to spontaneous emission by atoms excited by 
the detuned probe beam? The rate of this loss (ignoring 
reabsorption) is 7 x (excited population). We would like 
the ratio of this loss rate to the output loss rate Kfi to be 
small. In the optimal feedback regime and for A ^> 7, 
this ratio is given by 



7 M fi 2 /4A 2 _ AfxM 2AI sat 2AL 



KfJ 



- X' 



(4.15) 



For the typical values stated above, Eq. ( 4.15| ) is indeed 
small (~ 1CT 1 ). 

Another practical question is, how realistic is the zero- 
time-delay assumption for the feedback? It can be shown 
using the techniques of Ref. J52j that this assumption is 
justified providing the feedback delay time is much less 
than the lifetime of the trap due to the output cou- 
pling. If recent experiments ||3|, |9| are a useful guide, trap 



lifetimes are of order 10 s 



Feedback much faster 



than this should not be a problem. In fact, the time delay 
could be completely eliminated by feeding forward rather 
than feeding back. Linewidth reduction can be achieved 
equally well by controlling the phase of the atom field 
once it has left the trap as by controlling it inside the 
trap (but, of course, an integrated, rather than instanta- 
neous, current would be used for the control). 



monochromatic with small intensity and phase fluctu- 
ations. We used the normalized first-order coherence 
function g (1) (i) = G (1) (t)/ (7) Q as a measure of the 
phase fluctuations. As t increases, Ig*- 1 ^)! decreases from 
unity as the phase of the field becomes dccorrclatcd from 
its initial value. Its decay is characterized by the co- 
herence time r co h = i J °° I*/ 1 ) (t)\dt, or by its recipro- 
cal, the linewidth £. G^\t) also determines the output 
power spectrum, where the peak spectral height is given 
by 4(1) 1 1 regardless of lineshape, while for a Lorentzian 
the FWHM is also equal to I. 

We examined the linewidth as a function of the dimen- 
sionless atomic interaction strength, \ — 4/iC / k, where 
k is the output coupling rate and C is the atomic self- 
energy. There are four distinct regimes: the standard 
laser linewidth for \ <C 1, a quadratic dependence on \ 
for 1 <C X *C -y//!, a linear regime for ^//Z <C X *C fJ 2 , 
and finally another constant regime when \ ^ l^ 2 - The 
second and third regimes have Lorentzian and Gaussian 
output spectra respectively. The last regime is a conse- 
quence of quantum revivals of |<7^(t]L which are a di- 
rect effect of strong atomic collisions |3q] . This leads to a 
complicated structure in the power spectrum with many 
peaks contained in a Gaussian-like envelope. 

An important condition for atom-laser coherence is 
that the phase fluctuations be small in a particular sense. 
This is equivalent to requiring Bose degeneracy in the 
output, i.e., that the linewidth I be much less than the 
the output flux (I) — K/i. From the results presented 
here, this means that the laser output is only coherent for 
interaction strengths satisfying \ ~ M 3 ^ 2 ) i- e -j somewhere 
in the third linewidth regime. Therefore, collisions will be 
a problem for atom-laser coherence, especially for BECs 
formed by evaporative cooling where collisions are the 
dominant mechanism. On the other hand, if the atom- 
atom interactions are strong enough, the laser output will 
exhibit the interesting feature of quantum revivals. 

We also show, expanding upon Ref. po| , that this 
linewidth broadening can be significantly reduced by a 
QND feedback scheme. Basically, by feeding back the 
results of a QND measurement of the number fluctua- 
tions to control the condensate energy, it is possible to 
compensate for the linewidth caused by the frequency 
fluctuations. The very number-phase correlation created 
by the collisions is utilized to cancel their effect. We have 
shown that this linewidth reduction allows the output to 
remain coherent for interaction strengths up to \ — M 2 
rather than /i 3 ^ 2 , which is an improvement by a factor of 
y/Ji. For the reasonable parameters of C ~ 10 _2 s _1 and 
fi ~ 10 6 , this improvement in linewidth is of the order 
of 10 3 and in principle could increase coherent values of 
X = 4:(j,C/k from - 10 9 to - 10 12 . 



SUMMARY 
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APPENDIX A: SINGLE COMPLEX 
EXPONENTIAL APPROXIMATION FOR g (1) (t) 

For the majority of the calculations in this paper we 
assume that the first order coherence function can be 
given by 



g W (t) = \g^(t)\e 



(Al) 



That is, we assume that the complex nature of g^\t) 
(which is due only to collisions as shown in Sec. Ill D) 
is described by a single complex exponential for most 
regimes of the collisional nonlinearity. When the colli- 
sions are strong enough to cause revivals this approxima- 
tion necessarily breaks down. 

To determine the relative error caused by this approx- 
imation we continue the analysis of evolution due to col- 
lisions presented in Sec III D, i.e., for p = CcP- Specifi- 
cally, we quantify the relative difference between 



_ 1 

Texact — ~ / 6 




-/i(l-cos2Ci)^ 



(A2) 



and 



dt. (A3) 



This is done by multiple Taylor expansions [first of 
cos2Ct and exp(2iCi), and then expanding resultant ex- 
ponentials apart from exp(a£ 2 )] and using the equality 



t"< 



'dt = 



T[(n + l)/2] 
2{2pC 2 )( n + 1 )/ 2 ' 



(A4) 



where T[n] is the Gamma function. 

Now, T a pp rox differs from r GX act by both real and imag- 
inary terms. The dominant [for t < (/j,C 2 ) - 1/ ' 2 ] real term 
is 



(2pC-Qjf 
32 



1 



2 \pC 2 



3/2 



(A5) 



while the first three imaginary terms are 



(2/j.C - Q) 
8pC 2 



2 , 

r° 3 



8(pC 2 



. (A6) 

We can thus determine w by requiring that the imagi- 
nary terms vanish. This is essentially what the iterative 
procedure of Appendix B achieves. 

The series of imaginary terms above leads to a first 
choice of ujq = 2/iC, i.e., this sets the first term to zero 
and we are left with terms of 0[(^(C) _1 ]. To cancel the 
first two imaginary term we require 2/iC — Q = 2C/3, 
leaving terms of O^C)' 1 ]. Hence, w x = 2pC - 2C/3 
will largely ensure that T approx is real (since /j> 1). The 



question now is, how much different is r a p pr ox from r cxact 
using uj\l 

The dominant term in th e di fference is simply found 
by substituting u\ into Eq. (A5), giving 



1 



Mp 3 / 2 C V 2 



(A7) 



The relative size of this error term is obtained by com- 
paring to T cxact , which using the Taylor expansion equals 



1 r°° 99 1 

' e- 2 » ct dt = 



V/2C* V 2 



(A8) 



This corresponds to a relative error, given by (r approx — 
'?e X act)/'?c X act, of the order of /it . We have also verified 
the size of this error numerically for the full system dy- 
namics by simulations of |<7^(i)| and g^\t) exp(— iujit) 
(here uj\ is found in the same way as shown in Appendix 
B). These results also confirmed that the single complex 
exponential approximation breaks down when C ~ K7r/i 
(or x ~ 4-7T/X 2 ), as expected. 



APPENDIX B: ITERATIVE PROCEDURE FOR 
THE COHERENCE TIME r coh 

As stated in Sec. Ill B, the coherence time can be eval- 



uated numerically usi ng E q. (3.8), where the first guess 
for u> is given by Eq. (|3.9|). Subsequent corrections to Q 
are found by applying the procedure: 

T k = - Tr[a f (£ - iu k y 1 ap ss }/2 (a^a) , 
uj k+ i = oj k - Im(l/2Tfe). (Bl) 



At each step we have 
2r fc = 



g^(t)e- iUkt dt, 



(B2) 



which is simply a reexpression of Eq. (jBlJ) . Then using 
the assumption of Appendix A that 5 (1) (<) « \g^ (t)\e iQi , 
we have 



2r k 



\gM(t)\e*°- u *>dt. 



(B3) 



The simplest form for is a decaying exponen- 

tial, exp(— jt), which gives rise to a Lorentzian power 
spectrum as discussed at the end of Sec. II. In this case 



2r fe f 
and hence 



oo 



1 



e-^-^-^dt = -. 

o y-i{w-u>k) 



Im(l/2T fe ) « Lu k - lj. 



(B4) 



(B5) 



The more complicated form of |<7^(i)| = exp(— 7 2 f 2 ), 
which gives rise to a Gaussian power spectrum, also obeys 
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Eq. (B5), since in this case 



2r k 



-[-y 2 t 2 -i(uj-uj k )t]^ 



-(w- Wfe ) 2 /4 7 2 erfc 



2 7 

-l/i(u> - u k ). 



-i(u> ~ UJ k ) 



2; 



(B6) 



Thus, if the coherence function is in the regime where 
the output spectrum is Lorentzian or Gaussian, then 
u)k+i « w, and the first correction, will be sufficient 
for an accurate calculation. On the other hand, if the 
coherence function is in the regime of revivals then the 
approximation g^'{t) ~ \g^ 1 \t)\e lujt itself is no longer 
valid. Therefore, the above numerical method will be 
wildly inaccurate and we require another method for ob- 
taining the coherence time and linewidth. This is detailed 
in Sec. Ill D. 



Combining these leads to the following correspondence 
for the saturated gain term: 



d 



V[a']AW]- 1 p^-—n\ [ n+—n ] j Q(n,<p), (C3) 
which can be expanded in two different ways to give 



Equating these expressions leads to 



(C4) 



(C5) 



APPENDIX C: Q-FUNCTION 
CORRESPONDENCE FOR SATURATED GAIN 



l + }_d_Y 1 = + 

n dn ) \ dn , 



(C6) 



As stated in Sec. Ill C the master equation can be re- 
expressed as a Fokker-Planck equation for a convenient 
probability distribution. In this appendix we show that, 
for the Q function, the op erator correspondence for the 
gain term is given by Eq. (3.16). The individual super- 
operators in this expression have the corresponding dif- 
ferential operators, 



V[a)]p 
A[a)]p 



d_ 

On 



n Q(n, if), 



d 



n + T^n ) Q(n,ip). 



(CI) 
(C2) 



and hence 



- i } o(n '^ ) ' (c7) 



By then applying the Taylor expansion, (1 + x) 1 = 
J2kLo(~ x ) ■> we obtain the identity 

v^A^-'p -+ E(-£) Q( n ><p)- ( C8 ) 



fe=i 
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